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Abstract 

Recent observation of high energy neutrino in IceCube experiment suggests existence of superheavy dark 
matter beyond PeV. We identify the parent particles of neutrino as two degenerated right handed neutrinos, 
assuming the dark matter is the heaviest right handed neutrino. The 0{Vcb) ~ 0(10“^) flavor symmetry 
breaking accounts for the O(10~^) mass degeneracy of right handed neutrinos which is a sizable scale to 
explain the successful resonant leptogenesis at the PeV scale. At the same time, non-thermal production of 
the heaviest right handed neutrino gives the right amount of dark matter for Tuh ~ lOPeV. The footprint 
of flavor symmetry is left in degenerated mass spectra of extra Higgs multiplet and colored Higgs multiplet 
which may be testable for LHC or future colliders. 
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1 Introduction 

The standard model (SM) is a successful theory of gauge interactions; however, there are many unsolved issues 
such as how to generate Yukawa hierarchy, meaning of small neutrino mass, how to generate baryon asymmetry 
and what dark matter is. Recently, IceCube reported their observation of high energy neutrinos which might 
be a hint to solve some of those problems [1]. The most important problem of SM is how to stabilize large 
hierarchy between the electroweak scale, Mw ~ lO^GeV, and the Planck scale, Mp ~ 10^®GeV against quantum 
corrections. The elegant mechanism of stabilizing hierarchy is supersymmetry (SUSY) [2], which is the main 
target of Large Hadron Collider (LHC). The existence of light Higgs boson such as 125-126GeV [3] supports the 
idea of SUSY. 

In the minimal supersymmetric standard model (MSSM), as the Higgs superfields and are vector-like 
under the SM gauge symmetry Gsm = SU{3)c x SU{2)w x U{1)y, we can introduce ^-term; 

( 1 ) 

in superpotential. The natural size of parameter /i is 0{Mp), however /r must be 0{Mw) in order to allow 
electroweak gauge symmetry to break. This is so-called y^t-problem, which is solved by making Higgs superhelds 
chiral under a new U{l)x gauge symmetry. Such a model is achieved based on Ug-inspired extra U(l) model [4]. 
The new gauge symmetry replaces the /r-term by trilinear term; 

( 2 ) 

which is converted into effective /r-term when singlet S develops O(lTeV) vacuum expectation value (VEV) [5]. 
At the same time, the baryon and lepton number violating terms in MSSM are replaced by single G-interactions; 

GQQ + G^U^D^ + GU‘^E‘^ + G^^QL, (3) 

where G and G'^ are new colored superhelds which must be introduced to cancel gauge anomaly. Although 
these terms induce very fast proton decay, it can be suppressed by a S 4 havor symmetry [6] [7]. Therefore havor 
symmetry plays an important role in stabilizing protons in supersymmetric model. 

The iS '4 havor symmetry can also solve the hierarchy problem of Yukawa couplings and the problem of havor 
violating processes originated from the SUSY contribution simultaneously. The former can be achieved by not 
assigning the 5 ' 4 -triplet, on the other hand the latter can be done by assigning the 54 doublet for the hrst and 
second generations of left handed quarks. Therefore the hierarchies mu, rric <IC mt and m^, m^ <C mg are generated 
by the same manner as SU(2)w, as the mismatch of the sizes of representations among the helds in operator 
suppresses the coefficient of the operator. We should keep in mind that the hierarchy m(quark), m(lepton) ^ Mp 
is generated by the discrepancy of the sizes of representations of SU{2)w between left handed fermion and right 
handed fermion. This manner is also adopted for suppressing single G-interactions when G and are assigned 
to be 54 -triplets. From the sizes of the elements of Cabibbo-Kobayashi-Maskawa (CKM) matrix, the S 3 subgroup 
of 54 should be broken by 0(10“^). 

It is well known that heavy right handed neutrino (RHN) not only realizes small neutrino mass by see-saw 
mechanism but also generates baryon asymmetry of the universe through leptogenesis. If the first and second 
generations of RHNs are assigned to be 54 -doublet, then the mass degeneracy is solved to be achieved by 0(10“"^). 
Therefore successful resonant thermal-leptogenesis requires that the mass of the lightest RHN is O(PeV), which 
is in good agreement with the energy scale of high-energy neutrinos observed in IceGube. Furthermore, if we 
impose constraint on the reheating temperature, Trh < lO^GeV, for avoiding gravitino over production [8], the 
right amount of dark matter density is obtained through non-thermal production of the heaviest RHN. Although 
the flavor symmetry is broken, its remnant may be observed in the spectrum of extra particles such as G, G° or 
extra Higgs. 

2 Symmetry Breaking 

2.1 Gauge Symmetry 

We extend the gauge symmetry from Gsm to G 32111 = Gsm x U{1)x x U{1)z, and add new superhelds 
N'^, 5, G, G'^ which are embedded in 27 representation of Eq with quark, lepton superhelds Q, D^, L, E'^ 
and Higgs superhelds , E[^ . Where is RHN, 5 is Gsm singlet and G, G'^ are colored Higgs. The two U(l)s 
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Table 1: G 32111 assignment of superfields. Where the x, y and z are charges of U{l)x, U{1)y and U{l)z^ and Y 
is hypercharge. The charges of U{\),p and which are defined in Eq.(4) are also given. 


are linear combinations of U{1)^ and U{1)^ where Eq D S'O(IO) x U{\)^ D S'C/(5) x U{1)^ x 17(1)^, and their 
charges X and Z are given as follows 

X = Z =--Q^ + (4) 

The charge assignments of the superfields are given in Table 1. To break U{l)z^ we add new vector-like superfields 
which are originated in 351' -T 351' of Eq. The invariant superpotential under these symmetries is given 

by 


kU32lll 

= WoYWs + WgEW^, 

(5) 

fUo 

= -k Y^H°QD‘^ + Y^H^LE^ + Y^LN‘^ + 

(6) 

fUs 

= kSGG^ + XSH^H^, 

(7) 

Wg 

= Y'^'^GQQ + Y^^G'=U'=D‘^ + Y^^GU'^E'^ + Y^^G^QL + Y^^GD‘^N‘=, 

(8) 


= M<i,4>4>'=-k ^Y‘*’(4-4>n2, 

Mp ^ ' 

(9) 


where Mp = 2.4353 x 10^®GeV is reduced Planck scale and unimportant higher dimensional terms are omitted. 
Since the interactions Ws drive squared mass of S to be negative through renormalization group equations 
(RGBs), spontaneous U{l)x symmetry breaking is realized and U{l)x gauge boson Z' acquires the mass 

m{Z')^^gx{S), ( 10 ) 

where <C (5) is assumed based on the experimental constraint for Z'^ as follow 

, r 2.51TeV (ATLAS[9]) 

> I 2.26TeV (GMS[10]) ‘ ^ 

The constraint for Z' mass is not far from this bound. Note that the mass bound for Z' depends on sparticle 
mass spectrum [11]. In this paper we assume / {S) ~ 0(10“^). 

If = 0 in IT$ and the origin of the potential 4)“) is unstable, then <I>, develop large VEVs along 
the D-flat direction of (4>) = (4)“) = V, U{l)z is broken and U{l)z gauge boson Z" acquires the mass 


m{Z") = y^gzV. 

After the gauge symmetry breaking, since the R-parity symmetry defined by 

iir , 


R = exp 


20 


(fix — 8?/ -I- I5z) 


remains unbroken, the lightest SUSY particle (ESP) is stable. 


( 12 ) 


( 13 ) 
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Table 2: S 4 x Zg x x Zg x Zis x Z 2 assignment of superfields (Where the index i of the 5'4 doublets runs 
1 = 1,2, and the index a of the 5'4 triplets runs a = 1, 2, 3. The details of 54 are given in Ref. [13].) 


2.2 Flavor symmetry 


The superpotential defined in Eq.(5)-Eq.(9) has following problems. As the interaction Wg induces too fast 
proton decay, they must be strongly suppressed. The mass parameter M 4 , in 1T$ must be forbidden in order 
to allow U(l)z symmetry breaking. In Wq, the contributions to flavor changing processes from the extra Higgs 
bosons must be suppressed [12]. These problems should be solved by flavor symmetry. 

If we introduce 54 flavor symmetry and assign G, to be triplets, then Wg defined in Eq.(8) is forbidden. 
This is because any products of doublets and singlets of 54 do not contain triplets. Note that we assume full Eq 
symmetry does not realize at Planck scale, therefore there is no need to assign all superfields to the same flavor 
representations. In this model the generation number three is imprinted in G, G'^. Therefore they may be called 
”G Higgs” (generation number imprinted colored Higgs). 

As the existence of G Higgs which has life time longer than 0.1 second spoils the success of Big Ban nucle¬ 
osynthesis (BBN)[8], 54 symmetry must be broken. Therefore we assign to be triplet and to be doublet and 
singlet to forbid With this assignment, 54 symmetry is broken due to the VEV of $ and the effective 

trilinear terms are induced by non-renormalizable terms 


Wnrg = (GQQ + + GU^E^ + G^QL + GD'^N^). 

The size of effective coupling constants of these terms should be 


(^) m 

M2 


> 10 


-14 


to satisfy the BBN constraint [14]. The superpotential of gauge non-singlets $, is given by 

r4> 


V' 

, r^c\2 , l^c\2^ , ^ -‘2 r.T.2 , ^2^u,f.c\2 




= j^n[my + mr + mr]+ 


[<j>f + ‘i>2][wr + (‘j>2r + (‘&3)i 


2 \r^ 




Mp 


■[ 2 y 34 > i $ 2 ((< l >^)2 - ( ci >=) 2 ) + ($2 _ $ 2)(($-)2 + ($ 0)2 _ 2 ($ c ) 2 )] 


-k 


eYf 

Mp 


<i>3[y3$i(($^)2 - (d>“)2) + d>2((<l>=)2 + (<1>=)2 - 2($5)2)]. 


(14) 

(15) 


(16) 
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We assume the global minimum of the potential <i>°) is at S' 3 -symmetric vacuum as follow 


($1) = ($2) = 0, ($5) = (<1>^) = = M = (17) 

The assignments of the other superfields are determined based on following criterion. (l)The quark and 
lepton mass matrices reproduce observed mass hierarchies and CKM and Maki-Nakagawa-Sakata (MNS) matrices. 
(2)The third generation Higgs are specified as MSSM Higgs and the first and second generation Higgs 

superfields are almost inert Higgs (AIH). The representation of all superfields under the flavor symmetry 

is given in Table 2. 

In order to realize Yukawa hierarchies, we introduce gauge singlet flavon superfields Vi,XiY, Wi, Pw, Xp, P, X 
and fix the VEV of them by 



f 'm-susY 

Mp 

V Mp 

(W) 

2 

Mp 

e acw, 


1 

16 


= 0.1 = e, 


i^w) 

Mp 


{Xp) 

Mp 


C msusY ^. ^ 
\ Mp J 


{W 2 ) 

Mp 


e^Psw, 


Mp 


= 


(Pw) 

Mp 




0.01 = 

Mp 


msusY = O(TeV), 
(^2) 2 


(18) 


where cy = cos0y,sy = sin0y and Greek letters are complex except for e. In this paper, we assume original 
Lagrangian has CP symmetry and all parameters in it are real. Therefore the complex VEVs given in Eq.(18) 
induce spontaneous CP violation. We illustrate how to get these VEVs as follows. Without loss of generality, we 
can define (X ), {Xw) , {Xp) , (Vi) to be real by the redefinition of superfields. The superpotential of V is given 
by 


^ - ^2^) + ^(3ViV2 - V23)3. 

This superpotential and soft SUSY breaking terms give polynomial potential as 

V(s3) = Co + Cl S3 + C2S3 + C4S3, S3 = sin 30 y, 
where all coefficients are real. The minimum condition of this potential with 9 v, 

C3[Ci + 2C2S3 + 4 C 4 s|] = 0, 

has trivial solution C3 = 0 and non-trivial solutions which satisfy the equation 

Cl + 2C2S3 + 4 C 4 s| = 0. 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Which of these solutions is selected for the global minimum depends on the parameters in the potential. Since 
the solution C3 = 0 gives wrong prediction such as Vud = 0, we assume the solution of Eq.(22) corresponds to 
the global minimum. Whether the solution of Eq.(22) is real or complex also depends on the parameters. As 
we choose real solution, (V2) has no relative phase. On the other hand for V{Xp,P) and V{Xw,Wi,Pw), we 
assume {Wi) , (IU2), (P), {Pw) have relative phases which break CP symmetry. The scale of VEV is fixed as 
follows. If the superpotential of superfield 4' is given by 


W = 


nMp^ ’ 


(23) 


then the potential of 4' is given by 


U(4') = m||4'p 


A4'’" 

Mp^ 


|^pn-l|2 


(24) 


where ^ A ^ "msuSY is assumed. At the global minimum (4') ^ 0, as the each terms in the potential should 
be balanced, the scale of VEV is hxed by 


^ f msusY \ 

Adp Y Adp J 


(25) 
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2.3 Soft SUSY breaking terms 

We assume that the effect of SUSY breaking in hidden sector is mediated by gravity and induces soft SUSY 
breaking terms in observable sector. Since these terms are non-universal in general, large flavor-changing processes 
are induced by the sfermion exchange. From the experimental constraints on them, the assignment of quarks 
and leptons under the flavor symmetry are restrictive. 

After the flavor violation, soft breaking scalar squared mass matrices become non-diagonal. For the Higgs 
scalars, this gives the mixing mass terms as follow 

V = + m^e^ifffyiiy + mh^S;S2, (26) 

which compel extra scalars to develop VEVs as 

(ilf)=0(eV, ) = 0(e5)ud, (52)=0(e>., (27) 

where we put 

(H^) =vu = 150.7GeV, = Vd = 87.0GeV, {S 3 ) = u, > sjvl+v^ = 174.0GeV. (28) 

The VEVs of the extra scalar are very small, we call and ’’almost inert-Higgs” (AIH) and S 2 as ’’almost 
inert-singlet” (AIS). As Si does not develop VEV, we call it ’’inert-singlet” (IS). As the same effects affect the 
flavons, the VEV directions given in Eq.(17) and (18) are modified as follows, 

^ ($ 2 ) ^ 0(e4)U, ($-) = ^ + 0(e4)t/, ^ = eV + 0(e6), •••, (29) 

and so on. Note that dominant parts of scalar squared mass matrices of AIH and G Higgs are diagonal and 

degenerated. Due to the smallness of VEVs of AIH and AIS, the superpartners of AIH and G Higgs also have 
diagonal and degenerated mass matrices. Therefore the trace of Si symmetry is imprinted in their mass spectra 
which may be testable for LHG or future collider. 

2.4 The size of V 

The size of V given in Eq.(17) should satisfy following condition. The interactions of G Higgs are given by 

w + Q3Q3) = +G2 + G 3 )(U,^E^ + 21 / 3 ^ 3 ), (30) 

from which each G Higgs can decay to top and tau, where we assume G Higgs is lighter than G higgsino and G 

higgsino can decay to G Higgs. The decay width of G Higgs is given by 

r(G^f + r) + F(G^t + 5):^^ (5^+8x132), (31) 

where hig is mass of G Higgs and the factor 5 and 13 are the approximate value of renormalization factors [7]. 
As the the life time of G Higgs 


t(G) ~ 7.2 X 10 


-29 


TeV\ [ Mp 


sec, 


me j \ y 

should be shorter than 0.1 second in order not to spoil the success of BBN, we have to require 


me > 0.072 


V 


TeV. 


(32) 


(33) 


For 771(3 = ITeV, V/Mp > 0.52 x 10 should be satisfied. Assuming the A-term dominant potential as 

$^) = {-Ad>2[(d>(’)2 -h (4>^)2 -h (4>^)2] -h h.c.) + {D- term) + {F - term), (34) 
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the size of V is given by 


V 

Mp 


1 


—y 

Mp J 


X 10- 


- 7.5 


(35) 


which requires ~ 0.1 to give V/Mp = 10“®'®. Note that we assume Y = 0(1) means 0.1 < Y < 1 in this 
paper. Therefore the natural size of V is given by V/Mp < 10“®'®. In order to satisfy Eq.(33) for TeV scale G 
Higgs, the size of V should be in the region as follow 


0.52 X 10“®'® < — < 10“®'®. 
Alp 


(36) 


3 Quark Sector 

The superpotential of quark sector is given by 

W = H^QY^U^ + H^QY^D<^ + H^QY^ihU^ + (37) 

where Yukawa matrices are calculated as follows, 


J 

( Y^cye^ - 

-YYsye^ 

Y^yP^ \ 


Y^ = 

Y^sye^ 

Y^cye^ 

Yfsye^ 

(38) 

] 


0 

Ys^ J 


Y^ = 1 

( -Yi^C2ye® 

Yfs2ye^ 

Y^{asCy - ap^sy)e^ \ 


Y^S2ye^ 

Y^^C2ye^ 

Yf {asSy + ap-fCy)€^ , 

(39) 

\ 

K 0 

0 

Yfe^ J 



where C 2 v = cos20y and S 2 v = sin20y. As the Kahler potential receives the effect of flavor violation, superfields 
must be redefined as 


U‘^^Vk{U)U^, D^^Vk{D)D^, Q^Vk{Q)Q, 

/ 1 aie^ U2e'^ \ / 1 /3ie® /32e® \ 

VKiU)={ a*P 1 age® , Vk{D) = \ 1 /Sge® , 

V U2C^ a^e® 1 ) V /3|e® 1 ) 

( 1 qicve^ \ 

1 qisve^ , (40) 

gicye^ gisye^ 1 / 

in order to get canonical kinetic terms [15]. As the result, quark mass matrices are given by 

/ Yfycve® -Y^sve^ Y^cve"^ \ 

Mu = V^{Q)Y^VuVKiU) = Vu\ rfysye® Yj^cye^ Y^^sye^ , (41) 

V Y,^e^ alY,^e^ Yf ) 

( -n^C2ve® Y^S2vt^ Yi^aoe^ \ 

M'd = V^{Q)Y^VdVK{D)=Vd\ Y^Ds2ve^ Y^C2ve^ Y^ , (42) 

\ -qiCsvY^fP qiSsvY^^e^ Y^e^ ) 

ao = asCy - apjsv, Pd = a^sy + Qpjcy, jaup + |/3 dP = 1, (43) 


where some parameters are redefined for simplicity, such as Y^ + U 2 Y^ —>■ Y^ Note that the each elements in 
these matrices include only leading terms and the contributions to mass matrices from Yaih are negligible with 
this approximation. These matrices are diagonalized by the superfield redefinitions 


U^LuU, D^LdD, U^^RuU^^, 


RdD^ 


Cy Sy 

Ljr =1 —Sy Cy 

{Y^IY^)cy^ iJ^lY^)sy^ 


-{Y^jY^y 

0 

1 


(44) 

(45) 
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- 

Ijd — 


Ru = 


Rd = 


C2V 

S2V 


— S2V 
C2V 


~{Y4 lY^){aDC2V — PDS2v)i 


-{Yf/Yf)iaDS 2 V 


PdC 2 v)^ 


1 

Nua^e^ 

{Y3^IY3^)^^ ' 

-Nua3e^ 

1 

ose® 

{Y^/Y^)e^ 


1 

1 


{Q/Yi^y 


1 



-(6/>^3°)' 

A 1 


from which we get 

LlM[jRu 

lIm'jjRd 


diag(TO„, TOc, mt) = diag {{Y^ - Y^Y^ jY^)e^Vu, Y^e^Vu, Y^Vu) , 
diag(TOd, ms, mb) = diag {-Y^^e^Vd, Y^e'^Vd, Y^e^Vd) , 

Cay s^v Vub 


VcKM = l\jLd = —Say cay Vcb 

V Vtd Vts 1 

Vub = {rD{a*jjCv + /S^sy) - ru)e^ = {asTo - ru)(?, 

Vcb = tdUSdCv - al)Sv)e^ = aprDl*(?, 

Vtd = {rucsv — rD{aDC2V — /3i>S2y))e^ = [c3v{ru — asTn) + apTos^vl]^^, 
Vts = {ruszv — rD{aDS2V + PDC2v))e^ = [s3v{ru — agro) — apr£)Cay7]e^, 
tu = Y^IY^, ru=Y^lY^. 


The experimental values of CKM matrix and quark running masses at ITeV 


/ \Vud\ 

|K«| 

\Vub\ \ 

/ 0.974 

0.225 

0.355 X 10-2 \ 

\Vcd\ 

[Ucs] 

|U™| = 

= 0.225 

0.973 

4.14 X 10-2 

V \Vtd\ 


\Vtb\ J 

\ 0.886 X 10-2 

4.05 X 10-2 

1 J 


J = lm{VusVcbV:bV:s) = 3.06 x lO"® [16], 

m„ = (1.17 ±0.35) X 10-^ me = 0.543t[!;!]^^, m* = 148.1 ± 1.3, 

md = (2.401°;^?) X 10-3, m, = (4.9li;3) x 10-^ mb = 2Alt°olt (GeV) [17], 

are realized by putting the parameters at = Mp as follows 


ru — UsTd = 0.355, apTp = 4.14, arg(7) = 108.379°, 30y = 13.003°, 

\Y^ - Y^Y^/Y^\ = 1.5, 1 ^ 2 '^I = 0.71, 1^3^I = 0.28, 

\Y^\ = 0.38, [±2^1 = 0.78, |y3^| = 0.38, 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


(52) 


(53) 

(54) 


(55) 


where we use the renormalization factors given in [7]. As these parameters are consistent with the assumption 
that all factors Y are 0(1), quark mass hierarchy is realized without fine tuning. The soft SUSY breaking squared 
mass matrices of squarks are given by 





0(1) oie^ 
ale^ 0(1) 

4 * 5 

r\i. c'-‘ 


0 - 2 ^ 

0 ( 1 ) 



0(1) 

/3ie' 

/32e® \ 

Trii-\ i 

f 1 


Cye2 

/3re® 

0(1) 

/ 33 c® 


7rc' 

1 

Sye2 

/32*e® 


0(1) J 

\ 

V cye2 

Sye2 

0(1) 


and the squark A-term matrices are given by 


(56) 


V D 
Au = 

Ad = 


-VuUAuU’^ - VdDAoD^^ + h.c., 


/ A^Cyt^ — A^sye^ A^cye^ \ 

I Afsye® A^cye^ A^sye^ ) , 

V A^e^ alA^e^ A^ ) 


/ -Afc2ye3 Afs2ye"‘ Afa'e'^ \ 

I Afs2ye3 Afc2ye'^ A^P'e^ , 

\ -QiCayAfe^ giSayAfe® Afe^ J 


(57) 

(58) 

(59) 


where these matrices are defined for canonically normalized superfields. After the diagonalization of Yukawa 
matrices, the squared mass matrices are given by 


{'mu)scKM 


{'m\,)sCKM 


{'rnQ)sCKM 





0(1) 





0(1) 


(60) 

£4 


0(1) ) 


0(1) 

£5 


i 


0(1) 

£4 

(61) 


£4 

0(1) ) 

1 

/ 1 

£4 

\ 



1 


(62) 


£2 

0(1) / 



and the (1,1) elements of Ajj and Ajj remain real except for the next leading terms which are at most 0(e^°) 
and 0(e®), therefore the SUSY contributions to electric dipole moment of the neutron are negligible. Note that 
the (1,2) and (2,1) elements of {mQ)scKM are complex. The off-diagonal elements of squark mass matrices 
contribute to flavor and CP violation through the squark exchange, on which are imposed severe constraints. 
With the mass insertion approximation, the most stringent bound for the squark mass Mq is given by ex as 


/Im[(m|^ 


Q)i2irnjy)i2\ 


= € 


4.5 


< 4.4 X 10 


-4 


Mq 

TeV 


Mq > 72GeV, 


(63) 


where Mq = M(gluino) = M(squark) is assumed [18]. This bound is very weak and the SUSY FCNC problem is 
solved. Note that if CP symmetry was not imposed, Af* would be complex and the constraint on neutron EDM 
would gave stronger bound, Mq > 625GeV. 

The contribution to FGNG from AHI exchange is also suppressed due to the hierarchical structure of 
Assuming the mass degeneracy of GP-even AIH and GP-odd AIH, the strongest mass bound for AIH is given by 
D° - D° as niAiH > 79GeV [12]. 


4 Lepton Sector 


The superpotential of lepton sector is given by 

W = + H^LY^N'^ + LY^jhN^^ + 


where 


Y^ = 


yN ^ 


yM ^ 


Yfe^ 0 0 

Y^aEf^ —Y.^j3sw^^ Y^acwf^ 
Yfpse^ Yfacwe^ Y^ 

Y^^sv 


Y.^cv 


Yi^ + Yi^psw Yi^acw - Yl^-/w 0 


Yi^acw + F/7W Yf - Y,^Psw 

Y^fA’ 0 0 


0 


yMg6 


0 

rMA 


0 

0 0 Yfe^ 

OiE = bsCv -b bu,{/3cvSw + OiSyCw) — bpjwsv, 

Pe = bgSv + hyj{acvcw — Psysw) + bp^wcv, |« e |^ + \Pe\‘^ = 1 - 


(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


Because the all elements of MNS matrix are 0(1), the contributions of flavor violation in Kahler potential and 
Yaih and YM are negligible. Therefore the lepton mass matrices are given by 


/ Yfe^ 0 0 

Me = Y^vd = Vd YfaEe^ -YfPswe^ Y^acwP^ 
\ YfPEe^ Yfacwe^ Y^Pswe^ 


(69) 
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'M,, = 


Y^u 


Y.^e^Vu 

Xi = 


XiCy XiSv 0 

X 2 + X5/3sw x^acw — 2^47^ 0 
x^acw + Xijw X2 — X5/3sw 0 

(* = 1 , 2 , 4 , 5 ). 


y/WV 

The charged lepton mass matrix is diagonalized by superfield redefinitions 




Re = 


E —y LeE 


1 

^4 


1 

,2 



LeN 


ReE^, N 
1 0 0 
0 Psw —OLCW 

0 a* cw P* sw 


from which we get 

L^MeRe 


= diag(me, nir) = diag {Y^e^Vd ,+ a^c%^)e^Vd, Yfe'^Vd) , 


M' = 



-^31 + -^32 


Xi{M 21 CV + M22SV) Xl{M‘iiCV + M32SV) 
Xi{M 2 iCv + M22SV) 1^21 Y M22 M 2 iM^\ + M22M32 

xi{M^iCv + M32SV) M21M31 + M22M32 

M 21 = PswX 2 + (/3^Sw^ - Q;^c^) 2:5 - 074^041/2:4, 

M 31 = a*cwX 2 ,+{a* P + af3*)cwSwX5 + l3*JwSwX4, 

M 22 = —acwX2 + 2a/3c4i/S4i/a:5 — fi^wswXi, 

M32 = P*SWX 2 + (c^ - s^)2;5 - 0*744/044/2:4. 


The experimental values of charged lepton running masses at ITeV 

me = 4.895 X 10“^, = 0.1033, rrir = 1.757, 

are realized by putting the parameters at fi = Mp as follows 


(GeV) [17], 


\Y{^ \ = 0.30, 


|l^2 ^ (q:^c^ + /3^s^)| — 0.62, |l^3^| — 1-1, 


(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 

(81) 

(82) 


where we use the renormalization factors given in [7]. Charged lepton mass hierarchy is realized without fine 
tuning which is the same as for quark sector. 

In order to realize neutrino mass scale rrii, ~ O(0.01)eV, the relations 


2 _ {Yi^^'^vu? 

•La — 


YfV 


O(0.01)eV 


YfV - lO^^GeV, Ml = e^Y^^V - lO^GeV 


(83) 


are required. After the diagonalization of charged lepton Yukawa matrix, the squared mass matrix and A-term 
matrix are given by 


{'^\)sMNS — {Le'<tx\L*p) — ' 


0 ( 1 ) 

^4 


(84) 


Aie® O(e^) 0(e®) 

{^'e)smns = L^AeRe = Vd ( 0(e®) O(e^) O^e"^) 

0(e5) 0(e3) 0(e2) 


(85) 


where the (1,1) element of A'^ is real at leading order and the SUSY contribution to electric dipole moment of 
the electron is negligible. Based on consideration of the lepton flavor violations, the most stringent bound for 
slepton mass AIp is given by /i —?► e + 7 as 


—< 3 4 X 10"® ( \ 

Ml ■ ^OOGeVy 




Ml > SOOGeV, 


( 86 ) 
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where = M(slepton) = M(photino) is assumed [18]. Without CP symmetry, the constraint on electron EDM 
would give stronger bound, Ml > 1765GeV. 

The contribution to lepton flavor violation from AHI exchange is also suppressed due to the hierarchical 
structure of as same as quark sector. Assuming the mass degeneracy of CP-even AIH and CP-odd AIH, 
the strongest mass bound for AIH is given by ^ e -I- 7 as rriAiH > 38GeV [12]. 

For the canonically normalized superfields, RHN mass matrix is given by 


Mat 


V^{N)Y^VVKiN) = 

/ 1 0 

Vk^N) = 1 0 

\ 0 0 1 


Ml Me^ 0 \ 

Me^ Mi(l-fe4) 0 , 

0 0 Mg / 


M Ml, M3 = YfV€^, 


(87) 

( 88 ) 


which gives degenerated mass spectrum of RHNs as follow 


Ml ~ M 2 = Mi(l -b e^) < Mg 


Sn 


M 2 — Ml 
Wi 


(89) 


The right handed sneutrinos have same spectrum. In the early universe, the out-of-equilibrium decay of rig and 
Nl generates B-L asymmetry which is transferred to a baryon asymmetry by EW sphaleron processes. Following 
Ref. [19], the baryon asymmetry is given by 


Bf 


Keep 

3g* 


where g* = 340 is the degree of freedom of radiation, k is dilution factor which is given by 

1 


K = 


and ecp is given by 


a: In AT’ 

r(Mi) 


2 mM,y = = 


90M|, ’ 




1 Im(Ari2) f 2y/x ^ H-a;\ 

ecp = -x-^^ + v^ln- 

ZTT Kii \X — 1 ^ J 

'0.7PeV\ 




2'kKii5n‘ 


Ml 


From the order estimations as follows 


7^12 - K 


11 e , 


K : 


ecp 10“®, Ml ~ lO^GeV, 


(90) 

(91) 

(92) 

(93) 

(94) 


Ml J ’ 

we get the correct amount of baryon asymmetry Bf 10“^°. Right sign of baryon number corresponds to 

0° < 2BcP < 180°, Bcp = arg(ifi2). (95) 


5 Dark Matter 


The lightest SUSY particles of this model is singlino S2 which is the superpartner of AIS S'2. The superpotential 
of Higgs sector is given by 

W = AieS3(Hl^Hi^ + H^H 2 ^) + A3S3H^Hf 

+ A2e652[cy(Mf Mi" + H^H^) + sy(A7f H^)] 

+ Ue^evH^ + svH^)S 2 H^ + X^e^ S 2 [cv + syB^), (96) 


which gives 


(e^A4i;d)(e^A5)i;, 

eAiUs 


leV. 


M(s2) 


11 


(97) 















Although S 2 is not the dominant component of dark matter, it may help to explain the delay of structure 
formation [20]. The massless singlino si and LSP S 2 behave as extra neutrinos and change the effective neutrino 
generation number to 


3.194 [21], (98) 

where mz' < 4700GeV is assumed. This extra contributions soften the discrepancy of the expansion rate Hq 
between the measurements of type la supernovas and Cepheid variable and CMB data [20]. 

The interaction of bino is given by 


c = 

where higgsino has mixing mass term with S 2 as follow 

C = A 4 e® (Hj^) ~ e-^VuS2h^- 

Therefore bino life time is calculated as follow 


47r \msuSY J 


10-i°eV 


T ~ 10 ^sec, 


(99) 


( 100 ) 


( 101 ) 


which is consistent with standard cosmology. 

Five of the six flavon multiplets <f>, have lighter masses than lOOTeV and are produced non-thermally 
through the U{l)z gauge interaction. As the lightest flavon (LF) is quasistable, therefore it is the candidate of 
DM. Solving the Boltzmann equation with the boundary condition nLF{TRH) = 0, we get a relic abundance of 
the LF as [22] 


= 5.0 X 10-3 


Trh 

lO^GeV 


3 



The heaviest RHN n§ behaves like LF and has a relic abundance given by 

3 ,r\ 4 




lO^^GeVy/ Ms 
V / I lO^GeV 


( 102 ) 


(103) 


where M 3 is defined by Eq.(87) and we put M 3 ^ lO^GeV and ~ 0.1 based on following reason. As 
M 3 ^ rriLF is satisfied, the contribution of LF to Udm is negligible. can decay through the interaction in 
Kahler potential as 


K = [X*NsiN^yXwWiSySsY + (Pvv-contribution)]g 4 = + ^ 2 n 2 )siss, (104) 

and has a life time given by 

giOM^ 

r(n 3 ^ + Si + S 3 ) ~ ^ 10-'‘3eV ^ t ~ lO^^sec. (105) 

The daughter particles 2 decay to v -\- H and give neutrino flux at ~ M3/6 ~ PeV. From the latest results 
of IceCube [1], the mass and life time of DM are around PeV and 10^®sec [23]. If we identify rig as DM, Mg is 
fixed at lO^GeV from the IceGube results. Note that reheating temperature should be Trh < lO^GeV to avoid 
gravitino over production. After the three parameters in Eq.(103) are fixed as follows 


0.52 X 10-®-3 < -r^ < 10-3-3, Mg ~ lO^GeV, lO^GeV > Trh > Ms, (106) 

Mp 

the right amount of DM, ~ 0{1) — 0(10) is realized. This is an interesting prediction of our model. Eor 

the reheating temperature Trh ~ lO^GeV, the thermal mass of flavon $3 is estimated as follow 


rriTi^a) ~ (^Trh ~ lOGeV, 


(107) 
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which is small enough to avoid symmetry restoration. If inflaton decays through Planck suppressed interaction, 
for example, such as 


then reheating temperature is given by 

Trh ^ v'MpfW ~ ’ 

which requires inflaton mass TO(dt) ~ lO^^GeV for our model. 


(108) 


(109) 


5.1 Numerical analysis 

Finally we calculate detection probabilities of neutrino flavors emitted through dark 
of Ki 2 . The neutrino mass matrix is given by 

{Ml)exp = C/M^5diag(0,TO2,m3)C/|^jvs> 

( C12C13 S12C13 5136“*'^ 

— S12C23 — Ci2S23Sl3e*^ C12C23 — Si2S23Sl3e®^ S23C13 

S12S23 — Ci2C23Sl3e*^ —C12S23 — Sl2C23Sl3e*^ C23C13 


on the basis that charged lepton mass matrix is diagonalized. From the experimental values ; 

as follows [16], 

sin^ 012 = 

0.308 ±0.017, 


( 112 ) 

sin^ 023 = 

0 437+0 033 
u.‘±o ( _o,023’ 


(113) 

sin^ 013 = 

0 0 2 34+0 0020 
u.uzoa_Q QQ3g, 


(114) 

S/n = 

1 30+0-38 
i.O»_Q 27 , 


(115) 

Ato 2 i = 

(0.753 ±0.018) X 10-^eV^ 


(116) 

ATO 32 = 

(24.4 ±0.6) X 10-'‘eV^ 


(117) 

we fix the parameters as follows 




^*12 

= 33.709°, 


(118) 

023 

= 41.381°, 


(119) 

013 

= 8.799°, 


( 120 ) 

mi 

= OeV, 


( 121 ) 

m 2 

= 0.867756 X IQ-^eV, 


( 122 ) 

m 3 

= 5.015277 X lO-^eV, 


(123) 

6 

= 250.0°. 


(124) 

By the phase rotations of lepton doublets La, we define the diagonal elements of {Ml)exp to be real and non¬ 
negative and the real parts of (1,2) and (1,3) elements of {M[,)exp to be non-negative. The same procedure is 
performed for which is given in Eq.(76). The matching condition = {M[,)exp gives seven equations with 

nine free parameters as follows 

4 >, 0 a = arg(a), 9 b = arg(^). 

0c = Ow, Xl, X 2 , X 4 , 

X 5 . 

(125) 


matter decay and argument 
( 110 ) 

I diag(l,e*‘^,l), (111) 


Note that two of nine equations in matching condition are automatically solved due to the condition detM' = 
det(M')ea;p = 0. In order to solve the equations, two constraints should be added by hand. We fix {(j),X 4 ) to 
solve the equations numerically and calculate branching ratios and argument of Ki 2 . The former is given by 


Be 

Bp 


r(ni —>■ l^eH) + r(n2 — I'eH) 

Wh) 

r(ni —>■ VpH) + r(n2 ^pH) 

Wh) 


_fa_ 

xl + 2{xl + xl + xl) ’ 

X2 + x'^ + xl + 2X2X5SW cos 9 b — 2X4X5CW COs{ 6 a 

xf + 2 (x 2 + x1+ X 4 ) 


(126) 
— ,(127) 
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B. 


T{vH) 

where we put r(n 3 —?► ni + si + S 3 ) = r(n 3 —>■712 + 51 + S 3 ) by hand. The latter is given by 

9cp = arg(itri 2 ) = arg [(xfcvsv + 2 X 2 X 5 CW cos 9a — 2 xiX^sw cos(0b — 9c)) 

+ i{2x\cwsw sin(0o — 9^) — 2 x 2 X 4 sin0c)] , (130) 

where we assume the contributions from AIH are negligible for simplicity. Due to the neutrino oscillation, the 
neutrinos emitted by RHN-decay change the flavor, therefore the detection probabilities of neutrino at detector 
are given by 


r(ni —> PtH) + r(n2 —> VtH) X2 + x^ + X4 — 2X2X5SW cos 9 f, + 2X4X5CW cos( 0 a — 9 c) 

T^ivH) xf + 2 (x 2 + x^ + xj) 

= r(77i —> VeH) + r(7l2 —> VeH) + r(7li —> Pfj,H) + r(n2 —> P^H) 

+ r(7ll —> V-cH) + r(7l2 —>■ I'tH), 


,(128) 

(129) 


PiPe) 

= 

Be(|C/el|^ + |C/e; 

CO 

+ 




+ 

B^(\Uc1?\U^4? 

+ \Uc 2 ?\Uf, 2 ? 

+ \UcZ?\U^^\^) 



+ 


+ \Ue2?\Ur2? 

+ iC/eaPlffran, 

(131) 


= 

Be{\Uel?\U^,4? 

P\Ue2?\U^,2? 

+ \Ue^?\U^^\^) 



+ 

+ \U^2\^ +\U^^\^) 




+ 

Br(\U^4?\Url\^ 

+ \U^,2?\U^2? 


(132) 

P(Pr) 

= 

Be(|C/einC^rlP 

+ \Uc2?\Ur2? 

+ ic/e3nc/.3n 



+ 


+ \U^2\^\Ur2t 

■ + \UMUr^\'^) 



+ 

R,(|C/,ld+|C/,2d+|C/r3d), 


(133) 


where we assume the baseline is longer than PeVlml ~ 10 "‘pc. As the equation 30y = 13.003° has three 
solutions, 

9v = 4.334°, 124.334°, 244.334°, (134) 

we calculate 9cp,P{vc^^^t) for each cases. The results are given in Table 3. The dependences of detection 
probability P{Pe,fj.,r) on parameters are weak. Depending on (f), both sign of sin(20c'p) are possible. 


9v = 4.334° 


<(> 

X 4 

9a 

9b 

9c 

9w 

Xl 

X2 

Xb 

29cp 

P(Ve) 

P(^^^) 

P(yr) 

0 

1.5 

7.8 

118.7 

115.3 

31.5 

0.43244 

1.421 

1.262 

203.4 

0.231 

0.377 

0.391 

60 

1.5 

20.2 

178.5 

266.2 

17.0 

0.51667 

1.037 

0.844 

124.0 

0.238 

0.375 

0.387 


9v = 124.334° 


(!> 

X 4 

9a 

9b 

9c 

9w 

Xl 

X 2 

X 5 

29cp 

P(Pe) 


Pi’^r) 

0 

1.1 

25.5 

281.5 

276.6 

27.1 

0.43244 

1.399 

1.643 

138.4 

0.235 

0.377 

0.388 

60 

1.1 

21.5 

157.4 

118.5 

19.8 

0.51667 

1.221 

1.155 

239.2 

0.234 

0.375 

0.391 

120 

1.1 

30.3 

290.5 

169.5 

22.2 

0.61111 

0.696 

1.048 

38.4 

0.255 

0.371 

0.373 


9v = 244.334° 


<(> 

X4: 

9a 

9b 

9c 

9w 

Xl 

X 2 

X 5 

29cp 

P(^e) 

P(Vtl) 

P(yr) 

0 

1.1 

18.9 

282.6 

101.6 

153.2 

0.43244 

1.343 

1.685 

136.2 

0.233 

0.377 

0.389 

60 

1.1 

26.8 

177.4 

294.9 

165.8 

0.51667 

1.369 

0.975 

247.4 

0.234 

0.375 

0.391 

120 

1.1 

27.4 

315.6 

355.2 

150.4 

0.61111 

0.892 

0.887 

31.2 

0.256 

0.371 

0.372 


Table 3: The detection probabilities of neutrino flavors and argument of Ari 2 . Without loss of generality, we can 
define (xi, X 2 , X 4 , X 5 ) > 0 and 180° > (9a, 9w) > 0°. All angles are given in the unit of degree and Xi are given 
in the unit of 0.1-\/eV. 


6 Conclusion 

In this paper we consider S 4 flavor-symmetric extra U(l) model with taking account of high energy neutrino 
flux observed by IceCube and obtain following results. If we specify dark matter is the heaviest RHN, the mass 
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scale PeV is understood as follow. The structure of CKM matrix requires 0(10“^) flavor symmetry breaking 
and the symmetry should be non-abelian to suppress flavor changing processes induced by sfermion exchange. 
Due to this symmetry, two lighter RHNs form S' 4 -doublet and have same mass. As the 0(10“^) flavor symmetry 
breaking solves the mass degeneracy of S' 4 -doublet RHN by 0(10“^), the mass of RHN should be lOOTeV for 
a successful resonant leptogenesis. The 3-body decay of heaviest RHN with lOPeV mass generates PeV energy 
neutrino through the following decay of lighter RHN. Non-thermal production of the heaviest RHN gives right 
amount of dark matter. The information about leptogenesis may be extracted from neutrino flux. 

Although the flavor symmetry is broken, as the footprint of it is left in the degenerated mass spectra of almost 
inert-Higgs and G Higgs and their superpartners, the existence of flavor symmetry may be testable for LHC or 
future colliders. 
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